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Abstract 

There appear to be three challenges that any theory of dark matter must face: 
(i) why is Qdm of the same order as £1 Baryons ? (h) what are the near solar 
mass objects (~ 0.5Mq) observed by the MACHO microlensing project ? and 



[ (iii) understanding the shallow core density profile of the halos of dwarf as well 



as low surface brightness galaxies. The popular cold dark matter candidates, 
the SUSY LSP and the axion fail to meet these challenges. We argue that in 
the mirror model suggested recently to explain the neutrino anomalies, the 
. mirror baryons being 15-20 times heavier than familiar baryons, can play the 

role of the cold dark matter and provide reasonable explanation of all three 
Q-f above properties without extra assumptions. 
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■ I. INTRODUCTION 

Dark matter constitutes the bulk of the matter in the universe and a proper understand- 
ing of the nature of the new particle that plays this role has profound implications not only 
for cosmology but also for particle physics beyond the standard model It is therefore not 
surprising that one of the major areas of research in both particle physics and cosmology 
continues to be the physics of dark matter. 

Apart from the simple requirement that the right particle physics candidate must have 
properties that allow it to have the requisite relic density and mass to dominate the mass 
content of the universe, it should be required to provide a satisfactory resolution of three 
puzzles of dark matter physics: (i) why is it that the contribution of baryons to the mass 
density (Q) of the universe is almost of the same order as the contribution of the dark 
matter to it ? (ii) how does one understand the dark objects with mass ~ O.5M observed 
in the MACHO experiment @, which are supposed to constitute up to 20% of the mass || 
of the halo of the Milky way galaxy and presumably be connected to the dark constituent 
that contributes to Q ? and, finally (iii) what explains the density profile of dark matter 
in galactic halos; in particular, the evidence that the core densities of galactic halos remain 
constant as the radius goes to zero. 
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There are many particle physics candidates for dark constituent of the universe. Gen- 
erally speaking, the prime consideration that leads to such candidates is that they yield 
the right order of magnitude for the relic density and mass necessary to get the desired 
Q, dm pa 0.2 — 1. This is, of course, the minimal criterion. It requires that the annihilation 
cross section of the particles must be in a very specific range correlated with their mass. 
The most widely discussed candidates are the lightest supersymmetric particle (LSP) and 
the Peccei-Quinn particle, the axion. The first one is expected to have a mass in the range 
of 100 GeV whereas the mass of the second would be in the range of ~ 10~ 6 eV. The present 
consensus seems to be that value of flcDM is around 0.2 — 0.3, with Q\ making up the rest 
of the energy density of the universe at the moment. Compare these values with Qb — 0.05. 
The CDM contribution and the baryon contribution to Q are roughly of same order. On 
the other hand, the nucleon mass is very different from the masses of either the axion or the 
SUSY LSP. So to understand within the SUSY or the axion models why Q B ~ VL DM) one 
needs to work in a special range for the parameters of the theory. In either of these pictures, 
the MACHO observations must have a separate explanation. 

Furthermore, in recent years it has been emphasized that the LSP and the axion may also 
have difficulty in explaining the observed core density behaviour of dwarf speroidal galaxies 
which are known to be dark matter dominated. The point here is that both the axions and 
neutralinos, being collisionless and nonrelativistic, accumulate at the core of galactic halos, 
leading to a core density p(R) which goes like R~ 2 rather than a constant which seems to 
fit data better 0-0]. We will refer to this as the core density puzzle p|-p!0[ . 

Spergel and Steinhardt |§ have recently revived an old idea Jy], that dark matter may 
be strongly self interacting, to resolve this puzzle. They argue that, for the right range of the 
parameters of the particle, it may lead to a halo core which is much less dense and hence in 
better agreement with observations. Specifically they note that if the dark matter particle 
has mean free collision path of about a kpc to a Mpc, then the core on this scale cannot 
"keep on accumulating" dark matter particles, since these will now scatter and "diffuse out" . 
For typical dark matter particle densities of order of one particle per cm 3 , this requires a 
cross section for scattering of a ~ 10~ 21 — 10~ 24 cm 2 . Furthermore the properties of the dark 
matter particle must be such that it must not allow for dissipation of the thermal energy 
via emission of light particles; otherwise,the galactic core would cool and lead to an increase 
in core density. If these considerations stand the test of time, a theoretical challenge would 
be to look for alternative dark matter candidates (different from the popular ones described 
above) and the associated scenarios for physics beyond the standard model. 

A class of models known as mirror universe models have recently been discussed. These 
are motivated theoretically by string theory and experimentally by neutrino physics. They 
predict the existence of a mirror sector of the universe with matter and force content identical 
to the familiar sector (prior to symmetry breaking) |l~2V|l5l. Symmetry breaking might either 
keep the mirror symmetry exact or break it. This leads to two classes of mirror models: the 
symmetric mirror model, where all masses and forces in the two sectors remain the same 
after symmetry breaking [H| and the asymmetric mirror model [O] where the masses in the 
mirror sector are larger than those in the familiar sector. The mirror particles interact with 
the mirror photon and not the familiar photon so that they remain dark to our observations. 
Since the the lightest particles of the mirror sector (other than the neutrinos), the mirror 
proton and the mirror electron (like those in the familiar sector) are stable and will have 
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abundances similar to the familiar protons and electrons, the proton being heavier could 
certainly qualify as a dark matter candidate. We will show in the next section they can 
indeed play this role. 

It has been shown Hl6| , p!7| that, consistent with the cosmological constraints of the mirror 
universe theory, the mirror baryons have the desired relic density to play the role of dark 
matter of the universe. The additional neutrinos of the mirror sector are the sterile neutrinos 
that appear to be needed in order to have a simultaneous understanding of the three different 
observed neutrino oscillations i.e. solar, atmospheric and the LSND observations. In fact, 
one view of neutrino oscillation explanations of these phenomena fixes the ratio of familiar 
particle mass to mirror particle mass thereby narrowing down the freedom in mirror sector 
parameters. If indeed sterile neutrinos turn out to be required, mirror universe theory is 
one of the few models where they appear naturally with masses in the desired range. If we 
denote the mass ratio m p i/m p = (, then a value of ( ~ 10 — 30 is required to explain the 
neutrino puzzles. What is more interesting is that for the same parameter range required to 
solve the neutrino puzzles, mirror matter can also provide an explanation of the microlensing 
observations |16j- in particular why the observed MACHOs have a mass very near the solar 
mass and are still dark. 



II. ASYMMETRIC MIRROR MODEL IN BRIEF 



Let us start with a brief overview of the mirror matter models |P^,I1|. The basic idea 
of the model is extremely simple: duplicate the standard model or any extension of it in 
the gauge symmetric Lagrangian and allow for the possibility that symmetry breaking may be 
different in the two sectors. (See table below) There is an exact mirror symmetry connecting 
the Lagrangians (prior to symmetry breaking) describing physics in each sector. Clearly the 
W's, 7's etc in each sector differ from those in the other as do the quarks and leptons. 
When the symmetry breaking scale is different in the two sectors, we will call this the 
asymmetric mirror model 
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The QCD scale being an independent scale in the theory 
could be arbitrary. We will allow both the weak scale and the QCD scale of the mirror 
sector to differ from those of the familiar sector [TfJ and assume the same common ratio ( 
for both scales i.e. < H' > / < H >= A' /A = (. 



m, d, e, v e 



u' 



d'. 



W, Z, 7 , G 



w, z 1 , 7 ', a 

0', V' R ,... 



-Gravity- 



It is assumed that the two sectors in the universe are connected by only gravitational 
interactions. It was shown in JTBIJT^] that gravity induces nonrenormalizable operators that 
generate mixings between the familiar and the mirror neutrinos. This is one of the ingredients 
in the resolution of neutrino puzzles. To get an idea of how this works, note that the lepton 
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LHLH LHL'H' 




operators induced due to nonperturbative gravitational effects have the form w 
and L H M L H . After spontaneous breakdown (i.e. < H >= v and < H' >= v', we get for the 
mass matrix mixing the first generation neutrinos from each sector to have the form (in the 
basis Vf., ui): 

(1) 

where ( = v'/v defined above. To solve the solar neutrino puzzle via small angle MSW 
solution, we will choose A ~ 1 and £ ~ 10 — 30. This gives the sterile neutrino mass of 
order 10~ 3 eV or so, choosing Mpi ~ 10 18 GeV, which is also in the right range to explain 
the solar neutrino puzzle. Note that these are meant to indicate that the model leads to 
numbers in the right ball park. Emboldened by this result, we will consider the asymmetric 
version and look at its cosmological implications. 

In discussing cosmology, we first note that both sectors of the universe will evolve ac- 
cording to the rules of the usual big bang model except that the cosmic soups in the two 
sectors may have different temperature. In fact the constraints of big bang nucleosynthesis 
require that the post inflation reheat temperature in the mirror sector T' R be slightly lower 
than that in the familiar sector Tr (define (3 = T' r /Tr) so that the contribution of the light 
mirror particles such as z/, V etc. to nucleosynthesis is not too important. This is called 



asymmetric inflation and can be implemented in different ways ||18|| . We will see that, if we 
want the mirror nucleons to play the role of the dark matter, we will need a definite value of 
(3 depening on the choice of (. This in turn will help us to predict a value for the equivalent 
extra neutrinos at the BBN epoch (i.e. 5N V ). 

Before detailed discussion, let us first note the impact of the asymmetry on physical 
parameters and processes. First it implies that rrii —> C m i with i = n,p, e, W, Z. This has 



important implications for the nuclear and atomic physics of the mirror sector |13| , |23|| . For 
instance, the binding energy of mirror hydrogen is ( times larger so that the recombination 
in the mirror sector takes place much earlier than in the visible sector. With f3 = T' r /Tr 
as above, the mirror recombination occurs when the temperature of the familiar sector is 
C//3T r where T r is the recombination temperature in the familiar sector. The mirror sector 
recombination takes place before familiar sector recombination; this means that density 
inhomogeneities in the mirror sector begin to grow earlier and familiar matter can fall into 
them later as in typical cold dark matter scenarios. 



III. MIRROR NUCLEON AS DARK MATTER 

One can also compute the contribution of mirror baryons to the mass density of the 
universe as follows: 

^/3 3 C (2) 

Here we have assumed that baryon to photon ratio in the familiar and the mirror sectors are 
the same as would be expected since the dynamics are same in both sectors due to mirror 
symmetry. Eq. (2) implies that both the baryonic and the mirror baryon contribution to Q 
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are roughly of the same order, as observed. This provides a resolution of the first conceptual 
puzzle. Furthermore if we take fig ~ 0.05, then Qb' — 0.2 would require that (3 = (A/Q 1 ^. 
From this one can calculate the effective 5N V using the following formula: 

5N U = 3/3 4 + ^ 4 (^) 4/3 (3) 

where the last factor (11/4) 4 / 3 is due to the reheating of the mirror photon gas subsequent to 
mirror e + e~ annihilation. For £ = 20, this implies 5N V ~ 0.6 and it scales with £ as C~ 4//3 - 
Thus in principle the idea that mirror baryons are dark matter could be tested by more 
accurate measurements of primordial He 4 , Deuterium and Li 7 abundances which determine 
SN V . 

Clearly to satisfy the inflationary constraint of VLtot = 1, we need fl\ ~ 0.7. These kinds 
of numbers for cold dark matter density apparently emerge from current type I supernovae 
observations. It is interesting to note that if one were to require that Vie dm = 1, the 
mirror model would require that ( be much larger (more than 100) which would then create 
difficulties in understanding both the neutrino data and the microlensing anomalies. Thus 
mirror baryons seem to have just the right properties to be the cold dark matter of the 
universe. 



IV. EXPLAINING THE MICROLENSING ANOMALY 



Let us now turn to show how the mirror model accounts for the microlensing observations. 
We start with the four equations of stellar structure: 

dP/dr = -Gp(r)M(r)/r 2 (4) 



dM(r) jdr = 4tt r 2 p(r) (5) 
L(r)/47rr 2 = -(lQ/3)a SB (T 3 / pK)dT/dr (6) 

dL/dr = Aixr 2 e(r)p(r) (7) 

where n(r) is the opacity (cross section per unit mass) at radius r, c?sb the Stefan-Boltzmann 
constant, L(r) the luminosity at radius r, and e(r), the rate of energy generation per unit 
mass at radius r. We will need three terms in the equation of state (below) taken one or 
two at a time: 

P = (p/m)kT + (4<j S b/3c)T 4 + {h 2 /2m e )(3/87r) 2 (p/m) 5/3 (8) 

where the three terms represent gas pressure, radiation pressure, and (non-relativistic) 
degenerate electron pressure, m is the nucleon mass, m e that of the electron. We have 
neglected such niceties as keeping track of how many objects there are for each m of gas (2 
for H, 3/4 for He, etc) 
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We will make standard, illuminating if crude, approximations Jl9j in order to understand 
the ( behavior of the solutions to the above equations. First we write 



P = p GM/R, p = 3M/4tt R 3 (9) 

where P and p are roughly core averages. Here M and R are mass and radius of the star. This 
gives the useful relation 

P= (4tt/3) 1/3 GM 2 V /3 (10) 

To find the maximum mass of a (main sequence) mirror star, which is of interest to us, we 
note that as the mass of the star gets bigger, the core temperature rises. Therefore, of the 
three terms in the expression for the pressure in the Equation [8], we expect P g and P r to 
dominate. Following Phillips |19" |, we parameterize them as fractions of the total pressure 
P as below: 

P g =(3P, P r = (l-P)P (11) 

We eliminate T and solve for P from this parameterization to obtain 

PP = [(p k/m)\f3' 1 - l)/(4aW3c)] 1/3 (12) 

Using Equation fllCf ) again then gives 

M max ~ [(1 - f3)c/a S B] 1/2 G- 3 / 2 (k/m) 2 /f3 2 (13) 

As (3 approaches 1, the energy density is increasingly dominated by photons (relativistic 
particles) and stars become unstable. Taking a cutoff around (3 ~ 1/2 gives a maximum 
stellar mass around 70M Q . Thus the range for stars is roughly 0.07M© to 70M Q . From 
Equation ([H|) one sees, in the approximation that instability sets in at the same (3 indepen- 
dent of (, that M max varies as (~ 2 . For ( = 15, we get for the maximum mass of the mirror 
star O.5M , which is of the same order of mass as the MACHO microlensing events. Our 
model therefore provides a resolution of the microlensing anomaly that avoids the strong 



constraints of Freese et al |2(J for familiar sector white dwarfs. 

We want to point out here that we do not expect all of mirror dark matter to condense 
to form mirror stars. Instead, we would expect it to be in the form of a mixture of mirror 
dust and mirror stars. In this connection, it has been noted recently plj that current upper 



limits to scattering optical depths for Thomson scattering in early universe suggests that 
compact objects of any kind cannot be the main dark matter constituent. This would also 
suggest a mixed picture of the kind mentioned. 



V. SELF INTERACTION OF MIRROR MATTER AND HALO CORE DENSITY 

PROBLEM 

As noted in the introduction, there appear to be indications from the core density profile 
of dwarf and low surface brightness galaxies |5||| that the dark matter may need to be en- 
dowed with a significant self interaction. According to the analysis of Spergel and Stenhardt 
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||, the self interaction must be such that the collision cross section of dark matter particles 
should be of order 10 -23 (m/GeV) cm 2 for a one GeV particle corresponding to a DM mean 
free path of order of X ss = 300 kpc in a gas with number density given by Pdm/ 171 - Also we 
note that the cross section would scale linearly with the mass of the dark matter particle. 
The mean free path requirement cannot be met by the neutralino or the axion but is met 
quite naturally by the mirror dark matter forces |22| . 

Two obvious kinds of self interactions are the self interactions due to inter-atomic forces 
and nuclear forces. The latter are not very effective as can be seen by the following crude 
argument. From pure ( scaling, we can infer that un'n' ~ &nn(~ 2 - The value of C = 20 
would put the mirror nucleon cross section to be of order 10~ 26 cm 2 , which is much too 
small. Note that for ( = 20, one would need a cross section of order 10~ 22 cm 2 to get 
A = X ss . On the other hand, one would expect that there can be scattering of the dark 
matter particles (mirror atoms) due to inter-atomic forces. A crude estimate of such cross 
sections is given by cth>h> ~ tt {a' ) 2 , where a' is the Bohr radius of the mirror hydrogen 
atom. One can then estimate the mean free path of the dark matter particles in the mirror 
model to be Xdm ~ 3n( 3 x 10 16 cm. For ( = 20, this gives \dm ~ 0.3 kpc, which is not 
far below the required values for explaining the halo core density profile. Note that Spergel 
and Steinhardt || require Xdm ~ kpc to mpc. On the other hand, one could take the value 
of C = 30 and thereby get A ~ one kpc. 

We thus feel that mirror dark matter presents the best scenario for understanding the 
halo core density profile using a self interacting dark matter model. Of course more detailed 
numerical work is needed to confirm these qualitative conclusions. 



VI. STRUCTURE FORMATION 



Finally, we use the remainder of this article to make plausible that, in spite of the ( 2 
decrease in cross sections for most processes, the facts that (a) structure formation begins 
earlier for the mirror sector (because recombination occurs before matter-radiation equality) 
and (b) the higher mirror temperatures for the same processes, than familiar temperatures, 
permit formation of galactic and smaller structures. In doing this, we will make use of our 
previous work in [23| and [l6|, as well as that of Tegmark et al |f24 |. 



Much of the work of [2^] can be carried over to the present work, after suitable mod- 
ification to take into account the fact that, in the current model, the proton mass scales 
as (. Here, we will assume that primordial perturbations are "curvature" or "adiabatic" 
perturbations. This means that the scale of the largest structures are set by mirror Silk 
damping |2"5 |. 7' diffusion wipes out inhomogeneities until the 7' mean free path, 



A' = [a T C 



n 



(14) 



where (~ 2 <Jt is the mirror matter Thomson cross section and n' e is its electron number 
density, becomes greater than one third the horizon distance (ct). Silk damping turns off 
because the A' increases as z~ s while ct only increases as z~ 2 . 

First, we compute, from Silk damping, the masses of the largest structures in this picture. 
Structure formation starts with mirror sector particles, and familiar sector particles later 
fall into these. For numerical values below, we will take, h = 0.7 and Qb' = 0.2. We pick 
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t ~ (zi/z) 2 s with z\ = 4 x 10 9 and n B > = &b' p c z 3 / (C m p) with p c = l-9h 2 x 10 29 g/cm 3 . 
Silk damping stops at around z s d ~ 8( 3 which gives 

\ sd ~ 2.5 x 10 27 /C 6 cm (15) 
M sd ~ 10 54 /C 9 gm 

Note that, for ( ~ 10, this is about the mass (and size) of a large galaxy. This coincidence 
could be an important factor in understanding galaxy sizes should this model correspond to 
reality. 



As in ||23|| , we parametrize the separation of M s d from the expanding universe as taking 
place at 

Zstop = z M z sd (16) 

with 

Rc = Kd/z M (17) 

After violent relaxation we have for the proton temperature 

T p = GM G (m p /R G ~lO-*z M /( 2 ergs (18) 

with p, outside the central plateau, given by 

p(R) = A/R 2 ~ 10 26 z Af /(C 3 i? 2 ) gm/cm 3 (19) 

We now turn to the question of whether this isothermal sphere is likely to fragment and 
form mirror stars. For this we compute the amount of mirror molecular hydrogen since 
it is its collisional excitation (and subsequent radiation) that is believed to be the chief 
mechanism that provides cooling for formation of stars. If the rate for this mechanism is 
faster than the rate for free fall into the mass of the structure at issue, we can expect local 
regions to cool fast enough to result in fragmentation of that structure. We do here a very 
rough estimate of mirror galaxy fragmentation into mirror globular clusters, using the results 
of |p4| , but leave to a more detailed work further fragmentation into the O.5M structures 



predicted in fTE 



Reference |23] give a useful approximation to their numerical results for the fraction of 



molecular hydrogen, f 2 (fo denotes its primordial value): 

h(t) =f + (k m /ki)ln[l + x nht) (20) 

where, as a first try, k m can be taken as just the rate for H + e~ — ■> H + 7 (which they 
conveniently give as about 2 x 10 _18 T' 88 cm 3 s _1 ), while k\ is the rate for H + + e~ — > H + 7 
(2 x 10 _10 T _0 ' 64 cm 3 s _1 ). Equation (8) is the result of H 2 production from the catalytic 
reactions H + e~ — > H~ + 7 followed by H + H — > H2 + e~ competing against the 
recombination reaction that destroys the catalyst, free electrons, (approximately) as l/t 
(assuming constant density). Our goal here is to show from Equation (20) that it is plausible 
that f 2 , the fraction of molecular hydrogen, rises from its primordial value of 10~ 6 to the 
region above 10 -4 where cooling tends to be competitive with free fall. 
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First, we note that, if k =< vo >~ AT 1 cm 3 /s, for familiar e and p, we expect that, for 
mirror e and p, scaling with £ to go as ^~( 2+ ^> AT y , since a must go as (~ 2 , all factors of T 
must be divided by some combination of m e and M p , both of which go as ( (making this 
model much easier to compute from than that of ||23|| ). 

We now estimate fragmentation. From Equation (18) we see that the galactic temper- 
ature should begin at about 10 eV at a time when the cosmic temperature is about 1 eV 
and the cosmic gamma number density is about 10 9 / cm 3 . The rate for "compton cooling" 
is very fast at this high density (unlike at later times for the familiar case) and there should 
be rapid cooling to about 1 eV. We can now compute the Jeans length for fragments as a 
function of distance R from galaxy central. We use 

pj = (T/Gmf/M 2 (21) 

If we set the Jeans mass, M, to 4vrr 3 pj/3, we can solve for r obtaining (if we are careful 
to convert T in Equation (6) from ergs) 

r = R[10~ 7 ( 2 /z M ] 1/2 ~ 10~ 2 i? (22) 

Now inserting into Equation (20) gives the coefficient of the log term on the order of 
10 -35 and the coefficient of t in the argument varying from 10 -13 to 10 -17 as R varies from 
1 to 100 kiloparsecs while the free fall time ((Gp)^ 1 ^ 2 ) varies from 10 14 to 10 16 . This would 
appear to indicate the likelihood of fragmentation of the original Silk damping structure 
into smaller units, and the eventual formation of the 0.5M O( j oi black holes that explain the 



microlensing events of [16 



VII. CONCLUSION 

We have argued that the asymmetric mirror model |13[], originally proposed to solve 
neutrino puzzles and subsequently advocated |T6|JT7|] as providing an alternative dark matter 
candidate has the advantage of resolving the microlensing anomaly and possibly the core 
density problem of dark halos. Possible tests of these models are to narrow the allowed 
values of 5N V from more accurate observation of deuterium, He 4 and Li 7 and observing 
whether further accumulation of MACHO candidates lie in the mass between 0.1- 1 solar 
mass. Needless to say that if the underground searches for the cold dark matter now under 
way lead to a positive signal, mirror matter cannot be the dominant component of the dark 
matter of the universe. 

The work of R. N. M. is supported by the National Science Foundation grant under 
no. PHY-9802551 and the work of V. L. T. is supported by the DOE under grant no. 
DE-FG03-95ER40908. 
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